LOCALIZED POLYNOMIAL FRAMES ON THE BALL 



PENCHO PETRUSHEV AND YUAN XU 



Abstract. Almost exponentially localized polynomial kernels are constructed 
on the unit ball S'' in R"^ with weights W^(a;) = (1 - |x|2)M-i/2^ ^l > 0, 
by smoothing out the coefRcients of the corresponding orthogonal projectors. 
These kernels are utilized to the design of cubature formulae on B"^ with respect 
to Wi_i{x) and to the construction of polynomial tight frames in L'^ {B'^ ,W f_i) 
(called needlets) whose elements have nearly exponential localization. 



1. Introduction 

The construction of bases and frames on various domains, in particular on R*^ 
and on the d-dimensional cube, sphere, and ball, is important from many prospec- 
tives and has numerous applications. The example of Meyer's wavelets ^ and 
the ((9-transform of Frazier and Jawerth [6] clearly shows the advantage of using 
localized bases or frames for decomposition of function and distribution spaces on 
K'^ in contrast to other means such as atomic decompositions or Fourier series 
(in the periodic case). Three of their features, (i) infinite smoothness, (ii) almost 
exponential space localization, and (iii) infinitely vanishing moments, make them 
a universal tool for decomposing most of the classical spaces on R'', including Besov 
and Triebel-Lizorkin spaces. The key to this is that the coefficients in the wavelet 
or (/j-transform expansions precisely capture the information in the norms defining 
the corresponding spaces. 

Our primary goal in this article is to develop a similar tool for decomposition of 
weighted spaces of functions or distributions on the unit ball B'^ vaW^ {d > 1) with 
weights 

:= (1 - M>0, 

where \x\ is the Euclidean norm of a; G R*^. The situation here, however, is much 
more complicated than on M'' (the shift invariant case) or on the torus or even on the 
sphere due to several reasons: (i) there are no dilation or translation operators on 
B'^, (ii) the boundary of B'^ in combination with the weight W)j,{x) creates a great 
deal of inhomogeneity, (iii) orthogonal systems such as orthogonal polynomials on 
B'^ are much less friendly than the trigonometric system, and (iv) there are no 
uniformly distributed points on B'^ or on the d-dimensional unit sphere S"^. 

Our approach to the problem at hand will heavily rely on orthogonal polynomials 
in the weighted spaces L^{B'^ ,W^). The standard Hilbert space theory gives the 
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orthogonal decomposition 

oo 

(1-1) L\B'',W,) = Y,^Vt V^cUt 



n=0 



where Vf^ is the subspace of all polynomials of degree n which are orthogonal to 
lower degree polynomials in L'^{B'' ,Wf^). Note that dimV,f = ^ n''^^ so 

is a large subspace of L^. The orthogonal projector Proj„ : L^{B'^ ^W/^i) 
can be written as 



(Proj„/)(x) = / f{y)Pn{W^;x,y)W^{y)dy. 

Here Pn(W^; x, y) is its kernel and 5^ is the normalization constant of W^, namely, 
/b^ Wfj,{x)dx. It is crucial for our further development that the kernels 
PniyV^] X, y) have an explicit representation |19j in terms of Gegenbauer polynomi- 
als (see (|4.ip - (|4.2p below). Now, (jl.ip can be rewritten in the form 

oo 

/-^Proj„/, feL\B\W^). 

ri=0 

Denote by Sn the orthogonal projector of L'^{B'^, W^) onto J2u=o ® ^i!^^ i-^- ^nf 
ELoProjn/- Evidently, 

(1.2) {Snf){x) = b^, f f{y)Kn{W^;x,y)W^{y)dy 
with kernel 

n 

(1.3) ir„(VF^;x,y) :=^P,(W^^;x,2/). 
Consider now the kernel 

oo 

(1.4) Lit(x,y)=^s(^)p,(W^^;x,y), 

obtained by smoothing out the coefficients in the definition of the kernel Kn{W^; x, y) 
in (|1.3p by sampling a smooth function a. One of our main results in this article 
essentially states that if a G C°°[0, oo) is compactly supported, then L!^{x,y) has 
almost exponential (faster than any polynomial) rate of decay away from the main 
diagonal y — x in B^^ x B'*. To state this result more precisely, let us introduce the 
distance (see (|4.6p ) 

(1.5) d{x,y) := arccos { (x, y) + Vl- |.^lVl - on B'' 
and set 

W^,{n;x) (^1 - \x\^ + n'^^ ^ , xgB'^. 

Then (see ^ for any fc > there exists a constant Cfc > depending only on fc, d, 
/i, and a such that 

(1.6) mx, y)\ < Ck ^)^Jy^;^ii + ^d{x, y))~' 

The localized kernels provide a powerful tool for constructing cubature for- 
mulae on B'^ with weights Wfj,{x), /i > 0, that are exact for all polynomials of 
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degree n, i.e. in H^, and have positive coefficients of the right size. It is an im- 
portant feature of our cubature formulae (see ^J5|) that for all ^ > the knots are 
obtained by projecting onto B'^ sets of "almost equally" distributed points on the 
upper hemisphere S'^ in the knots are in fact almost equally distributed on 

B"^ with respect to the distance d{-, •) defined in (|1.5p . Currently very few families 
of cubature formulae with positive weights are known on B"^, among them is the 
family of the product type formulae [T6l[T^. However, the knots in these formulae 
are not almost equally distributed. 

Most importantly, the kernels enable us to construct localized polynomial 
frames in L^{B'^, W^) which is our primary goal in this article. Our construction is 
based on a semi-discrete Calderon type decomposition combined with our cubature 
formulae on the ball from Jj5] If we denote by 4" = our frame on B'^, where 

X = Wj^QXj is an index set consisting of the localization points (poles) of the frame 
elements, then we have the following representation of each / e L^{B'^, VF^): 

1/2 

The above clearly indicates that ^E* is a tight frame for L^{B'^ ,Wf^). The most 
important feature of the frame elements "0^ is their almost exponential localization: 
For C e Xj (the jth level in X) 

Ojd/2 

, -, Vfc > 0. 

Here the presence of the factor yJW^{2^;x) is critical; it reflects the expected 
influence of the boundary of i?'' and the weight W^{x) on the localization of the 
frame elements. Notice that the distance d{-,-) is also affected by the boundary 
of B"^. This localization of the ipi^ is the reason for calling them needlets. The 
superb localization of the needlets along with their semi-orthogonal structure and 
increasing (with the levels) number of vanishing moments enables one to utilize 
them for decomposition of spaces of functions or distributions on B'^ other than 
i^(i?'^, Wfj) such as L^{B''-, W^) (1 < p < oo) and the more general weighted Besov 
and Triebel-Lizorkin spaces on B''. We will report on results of this nature in 
a follow-up paper. 

These ideas were first used in |13j for the construction of frames on the unit 
sphere S"^ in R''+^. In [14] the spherical frames were utilized for decomposition of 
Besov and Triebel-Lizorkin spaces on the sphere. Further, this scheme has been 
applied in |15j for the development of frames on [—1,1] with Jacobi weights and 
then used in [8] for decomposition of weighted Besov and Triebel-Lizorkin spaces 
on the interval. 

This article is organized as follows. In ^J5]we outline the general principles which 
guide us in constructing localized kernels and frames on domains other than W^. In 
^we present some results on localized polynomial kernels on [—1, 1] with Jacobi 
weights. In prove our main results on localized polynomial kernels on B'^ with 

weights Wfj_{x), n > 0. In Sj5]we construct cubature formulae on B'^ with weights 
Wf_i{x). In f|6] we construct our needlet system and give some of its properties. 
Section [7] is an appendix, where we give the proofs of some results from the previous 
sections. 



(1.7) M^) 
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Throughout this paper positive constants are denoted by c, ci , . . . and they may 
vary at every occurrence. As usual the constants may depend on some parameters, 
which are indicated exphcitly in some important cases. The notation A ^ B means 
ciA< B < C-2.A. 

2. General principles for constructing localized kernels and frames 

Let (E^ /i) is a measure space with E a metric space and suppose that there is 
an orthogonal decomposition of L'^{E, 

oo 

(2.1) i2(i?,M) = E0^"' 

71=0 

where V„ is a subspace of dimension dimV„ ^ n^, 7 > 0. Let P„ be the kernel of 
the orthogonal projector Proj„ : L'^{E, fi) — > V„, i.e. 

(Proj„/)(x) = / P„{x,y)f{y)d^l, f G L\E,^l). 
Je 

Notice that Pn can be written in the form Pn{x,y) = Pj{x)pj{y), where 

{pj} is an orthonormal basis for V„. Then Kn X]j=o kernel of the or- 

thogonal projector onto X]"=o ® "'^f • most cases of interest the kernel Kn{x, y) 
has poor localization, examples include the trigonometric system, orthogonal poly- 
nomials in one or several variables on various domains. 

Localization principle. It seems to us that there is a general localization princi- 
ple, which says that for all "natural" orthogonal systems, if the coefficients of the 
kernel Kn are smoothed out as in (|1.4p by sampling a C°° function, then the result- 
ing kernel has "excellent" localization around the main diagonal y — x in E x E. 
To be more specific, suppose a g C°°(M), a is even, and a is compactly supported 
or a e 5 (the Schwartz class of rapidly decreasing C°° functions on R) . Define 

C30 

(2.2) L„(x,y) 

Then for all "natural" orthogonal systems, the kernel Ln{x,y) decays away from 
the main diagonal y = x at nearly exponential (faster than any polynomial) rate 
with respect to the distance in E. 

In the case of the trigonometric system this principle is well-known and widely 
used. It is a fundamental fact in Harmonic Analysis that the Fourier transform of 
every function / in the Schwartz space S = S{M.'^) belongs to the same space. Con- 
sequently, any trigonometric polynomial Ln{t) :~ X]"=-n ^i'^*'^* with coefficients 
{a^} obtained by sampling a compactly supported C°° function has faster than any 
polynomial rate of decay away from zero. To make this more precise, let 

where a is compactly supported and a € C°°(M). Then L„ is a trigonometric 
polynomial of degree an and one easily shows the following localization of L„: For 
any k > Q and r > there exists a constant > depending only on fc, r, and a 
such that 

(2.3) \L^:\t)\<Ckj^^^^. teh^,7r]. 
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This estimate will serve as a prototype for our further localization results. 

For Gegenbauer polynomials and spherical harmonics the localization principle 
is established and used in T3] and also follows by the general result in [7j on the 
spectral properties of elliptic operators. For Jacobi polynomials it is proved in 
[3] and [15] (see Theorem 13.11 below). For Herniite and Laguerre polynomials the 
localization principle is established in . We will establish it here for multivariate 
orthogonal polynomials in L^(i?'',W^) (see Theorem 14. 2|) . We believe that the 
localization principle is valid in more general settings as well. 

For our purposes we restrict our attention to "smoothing functions" a satisfying: 

Definition 2.1. A function a is said to be admissible ifaCz C°°[0,oo), a(t) > 0, 
and a satisfies one of the following two conditions: 

(a) suppa C [0, 2], a{t) — 1 on [0, 1], and < a{t) < 1 on [1, 2]; or 

(b) suppSc [1/2,2]. 



There are two important applications of the localized kernels L„(x, y) from (|2.2I 
(i) If a is admissible of type (a), then the operator 

{Cnf){x) := / L„{x,y)f{y)dfi{y) 



apparently satisfies: £„/ = / for all / e I]"=o®"^!^ and £„/ e Z^i^lo®"^"- 
These along with the superb localization of L„ (to be established) makes £„ a useful 
tool. We will see this operator at work in the construction of cubature formulae on 
the ball in %\ 

(ii) Kernels Ln{x,y) with a admissible of type (b) are a valuable tool for con- 
structing localized frames. Let, in addition, a satisfy the conditions: a(t) > 
and 

(2.4) a^{t) +a^{2t) = 1, [1/2,1]. 

Then 



(2.5) ^0^(2-^0 = 1, te[l,oo). 

It is easy to construct such functions (see e.g. jlF). Define 

oo 

(2.6) Lo(a;,y) := Po(a;,?/) and Lj(x,?/) := ^a(^^^)p^.(x,?/), j = l,2,. 
and denote briefly 

{Lj*f){x):^ / Lj{x,y)f{y)dfi{y). 



One easily obtains the following semi-discrete Calderon type decomposition (see 

e.g. m) 

oo 

(2.7) f = Y,Lj*L^*f tor f eL2{E,fi). 

To get a completely discretized decomposition of L2{E,ii) one can use quadrature 
(cubature) formulae, if available. Assume that there is a quadrature formula 

(2.8) / /rfM ^ E ^«/(^) 
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with Xj C E and > 0, which is exact for aU functions / of the form f — gh with 

After these preparations we now define the frame elements by 

(2.9) t^^ix):=^^-L,ix,0 for ^ ^ X„ j = 0,1, . . . . 

The V^'s inherit the locahzation of the kernels Lj, which is almost exponential in 
all cases of interest. This is the reason for calling them needlets. 

We write X :— U°^QXj, where any two points (^,uj E X (from levels Xj ^ Xk) 
are considered to be different elements of X even if they coincide. We use X as an 
index set in the definition of the needlet system 

The next statement shows that is a tight frame in L^{E,fj,). 
Proposition 2.2. If f e L'^{E,fj.), then 

00 

(2.10) f^Y.Y.{f,^^)^^^J2{f,i;^)^^ mL^{E,fi) 

3=0 iex^ iex 

and 

(2.11) \\f\\mE,,) = {J2\{f,m'y^'- 

The proof of this proposition is a mere repetition of the proof of Theorem 4 in 
[T5] (see also [13]) and will be omitted. 

3. Localized polynomial kernels on [—1, 1] 

The Jacobi polynomials {Pn''^''}^^o constitute an orthogonal basis for the weighted 
space L'^{[-l,l],Wa,i3) with Wa,p{t) := {I - t)°' {I + , a,f3 > -1. We let Ca,i3 
denote the normalization constant of Wa,f3, i.e. :— J^^Waj3(t)dt. It is well 

known that [T7j 

where 

^ r(a + /3 + 2) r(n + Q + l)r(n + /3 + l) 

" ' " r(a + l)r(/3 + 1) (2n + a + /? + l)T{n + l)r(n + a + /? + 1) ' 
For / £ i^([— 1, l],Wa^ff) the Fourier expansion of / in Jacobi polynomials is 

fit) = ^ d„(/)(/.("-«)-iF("-«(0, dM) = Cc.,P / f{t)p.i''^^'Ht)wc.At)dt. 

The nth partial sum of this expansion can be written as 

{S,J)ix) - f]d,(/)(/i(."'«)-ip,("'«(a;) = c„,^ r f{t)KtP\x,t)w^.0{t)dt, 
3=0 -^-1 

where the kernel is given by 

(3.1) K(-^)(.,t)=^(/.;-^))"F,(-«(.)F,(-'''(,). 
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The grand question here is: What is the locahzation around the main diagonal 
?; = a; in [— 1, 1]^ of a polynomial kernel of the form 



oo . , 

(3.2) /.-^(.,y)=E<-)('^rT^^"'^(-W"''(^)' 



3=0 

where a G C°°l 

To address this question, denote 

Wc,,p[n-x) (l-a; + n-2)"+i/2(i + 2. + „-2),3+i/2^ 

Theorem 3.1. |15| Let a,(3 > —1/2 and let a be admissible according to Defini- 
tion \2.1\ Then for every fc > there is a constant Ck > depending only on k, a, 
[3, and a such that for < 0, (/> < tt 

(3.3) |L,^'^(cos6',cos0)| < cfe- 



\/'Wa,i3{n; cos9)^Jwa^i3{n] cos </))(! + n\9 — 
Here the dependence of Ck on a is of the form ct — c(a, /3, k) maxo<jy<fe ||a^''''||L=° • 



For the proof of this theorem it is important to establish estimate p.3p first in 
the particular case when </> = (the localization of L"'''(x, 1)). Set 

oo . , 

(3.4) L-^(x) := i-^(x, 1) = Y.a(l) (hf^^'Y P,(-^)(l)P,(-«(x). 

Since J7, (4.1.1), p. 58] 

p(«./3)n ^ ^ + = r(n + a + l) 

" ^ ^ V " / r(a + i)r(n+i)' 

it is easy to verify that 

where := r(/3 + l)/r(a + /3 + 2). 

Now the key role is played by the following theorem, which will also be critical 
for the proof of our main localization result (Theorem 14. 2p . 

Theorem 3.2. [31 [TS] Let a be admissible and assume that a > (3 > —1/2. Then 
for every fc > and r > there exists a constant Cfc > depending only on k, r, a, 
P, and a such that 



(3.6) ^LZ-P{cos0) 



<^kj-—^, 0<e<7r. 
(1 + nff)'' 



dx'' 

The dependence of Ck on a is of the form Ck — c{a, (3, fc, r) maxo<i/<fc ||a''^^||L°° • 

This theorem is proved in [3j with a admissible of type (a) and in (15j with 
a admissible of type (b). The proof in [TS] rests on the localization properties 
of trigonometric polynomials given in (|2.3p . while the proof in [3] is based on a 
property of Jacobi polynomials; it can be carried out with a admissible of type (b) 
as well. Estimate p.6p was proved earlier in [13] in the case a — P — X~ 1/2 (with 
A a half integer) and utilized for the construction of frames on the n dimensional 
sphere. For the reader's convenience we give the proof of Theorem 13.21 (following 
the idea from [3]) in the appendix. 
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Theorem 13. II is established in [l5j. Its proof rests on Theorem 13. 2 



4. Localized polynomial kernels on the unit ball 

It is known (see [19^) that the orthogonal projector Proj„ : L^{B'^,Wfj^) 
can be written as 



(Proj„/)(x) = / f{y)PjW^;x,y)W^{y)dy, 
where if /i > the kernel P„(W^;a;, has the following explicit representation: 
(4.1) P^^{W^-x,y) = hr'^^ 



C^: [{x, y) + u,/l^,/l^) (1 - u'T-'du, 



where (a;, y) is the usual Euclidean inner product, C,^ is the nth degree Gegenbauer 
polynomial, and 

d-l 

The case /i = is a limit case and we have 

\ + n 



(4.2 



-Pn(Wo;x,y) 



2A 



c,^((x,y) + yr^yT~H^) 

For an admissible a (according to Definition 12. ip we define 

oo 



J=0 



The explicit representation (|4.ip gives 



(4.3) = 
where is defined by 



i \ J + A 



n/ A 



cHx). 



Since 



, _ r(A + l/2) r(n + 2A) 

"^^ r(2A) r(n + A + l/2) " ^ 



it readily follows from p.Sp that i^^^^^'''' ^^^(x) = L^{x). Then by Theorem 
we get the following estimate: For all fc, A > and r > there exists a constant 
Cfc > depending only on fc, r, A, and a, such that 



(4.4) 



< Cfc 



-,2A+2r+l 



(1 + n6 



0<e <Tr. 



Distance on B'^. In order to show that Ll^ is a well localized kernel and for 
our further development, we need to introduce an appropriate distance in B'^ that 
takes into account the fact that B'^ has a boundary. In [18] it is shown that the 
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orthogonal polynomials on the unit ball and those on the unit sphere are closely 
related by the simple map 



(4.5) xeB'^^x' -.^ (x, vT~R2) e 5"^, 

which "lifts" the points from B"^ to the upper hemisphere S'l in that is, 

S'l := {x G S''' : Xd+i > 0}. This relation leads us to the following distance on B"^, 
which will play a vital role in the following: 

(4.6) d{x, y) := arccos | {x, y) + ^1 - k| Vl " l^p} ■ 



In fact this is the geodesic distance between x' := (x, a/I — |a;p) and y' (y, ^1 — |yp) 
on 5^ C R''+^ and, consequently, it is a true distance on B'^. This distance has 
been used to prove various polynomial inequalities, see the discussions in [2] and 
the references therein. 

The map (|4.5p also leads to a close relation between the spaces Lp{B'^, Wq) and 
LP{S'^, duj), where doj is the surface measure on S"^. This allows us to derive results 
on LP{B'^,Wo) from those on U'{S^^,dijj), which are also easier to prove. For these 
reasons we will prove our results only in the case /j, > 0. 

The following lemma provides an important relation between c?(-, •) and the Eu- 
clidean norm | • | in B'^ . 

Lemma 4.1. For x,y G B'^ , we have 

(4.7) I \x\ -\y\\< -^d(x, y) (VT^RF + 
and hence 

(4.8) I ^/l^W - < ^d{x, y). 

Proof. Let < a,/? < 7r/2 be defined from |x| = cos a and \y\ — cos/3. Using 
spherical-polar coordinates x = \x\$^ and y — \y\C,, where C G 5"'^^, we see that 

d{x,y) = arccos (cos a cos /3(^,C) -I- sin a sin /3) > arccos(cos(Q! — /?)) 

which yields d{x,y) > \a — /3\. On the other hand, since < a,/3 < 7r/2, we have 
cos > cos(7r/4) = \/2/2 and, consequently, 

sm a + sm p = 2 sm cos > v 2 sm . 

2 2 - 2 

Using the above we obtain 

II I I II I fli o • ■ ^ + P 

\\^\ ~ \v\ \ — \ cos a — cos p| = 2 sm — ^ — sm — ^ — 

1 , . . 1 



< ^|a-/3|(sina-|-sin/3) < —d{x,y){^l-\x\' + ^l-\y\^). 

Thus (|47l) is established. Estimate ([48ll follows immediately from ([iJ]) . □ 
Let us define 

(4.9) W,,(n;a;) := (yT~|^ + n-i)^^, x e B''. 

Our next theorem shows that the kernels L^' are almost exponentially localized 
around the main diagonal y = x va B"^ x B"^. 
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Theorem 4.2. Let a be admissible. Then for any fc > there exists a constant 
Ck > depending only on k, d, /i, and a such that 

y^Wf,{n;x)y^Wf,{n;y){l + nd{x,y)f 



Remark 4.3. TheoremlJTE as well as Theorems \3.1aKM can obviously be modified 
for the case when a £ {k sufficiently large) in place ofaG C°° . 

We will derive Theorem 14.21 when /i > from estimate (|4.4p and the foUowing 
lemma, using representation (|4.3p of i^. The proof in the case = is easier and 
wiU be omitted; it utilizes (|4.2p . 

Let us denote briefly 



(4.11) tix, y- u) := {x, y) + I - \x\^ ^ I ~ \y\^ . 

Lemma 4.4. Let 7 > —1, fc > 87 + 4, and n > 1. Then for x,y E B'^ 
(1 - u'^)'^du 



(4.12) 



(1 + n^l - t{x,y]u)Y 



< 



(v/1- |a;|2 + n-i)'r+i(^l-|2/|2 + n-i)^+i(l + nd{x, 2/))'^-3-r-4 ' 
where c > depends only on 7, k, and d. 
Proof. Denote briefly t :— t{x, y; u). Then we can write 



l-t=l-{x,y)- v/T^^v/T^7F + (1 - u)^/l^^l^, 
which implies 



(4.13) l-t>l-{x,y)-^l-\x\^^l-\y\^ 

= 1 - cos d(x, y)^2 sin^ ^^^EjII > .^d(x, y f 

and 

(4.14) l-t> \d{x, yf + (1 - u)^l - |x| Vl - \y? 



By (HT^ . we have 



>(i-^.)yr~N^Vi^ 



1 {\~u^)'^du 



^^■^^^ 7-1 (1 + nVT~tY - (1 + nd(x, 2/))'= 

Inequality (|4.12[) will follow from this and the estimate: 

(1 - u^ydu 
-1 (1 + n^JT^tf 

< 



(4.16) 



To establish this last estimate, we split the integral over [—1,1] into two integrals: 
one over [—1,0] and the other over [0, 1]. For the integral over [—1,0] we write the 
factor {l + n^/T^tf as the product of (l + nVT^)''"^''"^ and {l + ny/T~tf^+'^. 
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Then we apply inequalities (I4.13P and ()4.14p to the first and the second terms, 
respectively. This gives 



„ fO 

-1 

< 



- (l + nd(a;,y))fc-27-2 [n^ ^1 - \x\^ ^1 ~ \y\^{l - u)p+^ " 



We now estimate the integral over [0, 1]. Denote briefly A :— y/l^^Jx^ ^yl^^\y^ . 
Using (|4.14p and applying the substitution s = An'^{l — u), we get 

< c / , ^ ' du 

Jo {1 + ny'd{x,y)^ + A{1 - u))'^ 

^ c r^"' £2 

- (An2)7+i 7^ (1 + ^n2d(a;, yFT7)fe ^ 

< 



+ y))fc-27-3 _^ ^n2d(a;, y)2 + s)27+3 

„„-27-2 



< 



yl7+i(i + nd(a;,y))'^-27-3- 

Putting these estimates together gives ()4.16p . 

To complete the proof of Ij4.12p we need the following simple inequality (see 
inequality (2.21) in [15]): 

(4.17) {a + n-^){b + n-^) <3{ab + n-^){l + n\a-b\), a,b>0,n>l. 
Inequalities ([48]) and (|4.17p yield 



(yr^N^ + n-i)(v/T^^ + n-i) 



< 3(v/r^^v/l^^ + + nd{x, y)). 

This along with (|4T5)) and (|4T6| implies ((4T2)) . □ 



Proo/ o/ Theorem\4^ For t = cos 6*, < 6* < tt, we have 6'/2 - sin6'/2 - Vl^I. 
Therefore, estimate ()4.4p with r = is equivalent to 



^2A+l 

|L^(i)|<Cfe , -l<i<l. 

Now, (I4.10p follows readily by Lemma [4.41 □ 
The estimate of \LI^{x,y)\ from Theorem 14.21 enables us to control its LP-norm. 



Proposition 4.5. For < p < oo, we have 

(4.18) (/ mx,y)\m,{y)dy) < c(— -) 



X e B'\ 



Proof. If < p < oo this proposition is an immediate consequence of Theorem 
and Lemma [4.61 below, taking into account that estimate (|4.10p holds for an arbi- 
trary k. In the case p = oo estimate (14.181) follows by (|4.10p and (|4.8p (see also 
estimate (|422)) below). □ 
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Lemma 4.6. If a > d/p + 2n\l/p- l/2\, ^ > 0, < p < oo, then 

(4-19) Jp--= [ ..w ^^ < cn-^W^n; 

^ ' is. W^(n;y)f/2(l+nd(a;,y))^f - ' ' 

where c > depends only on p, /i, and d. 

Proof. Let > (the case = is easier). Three cases present themselves here. 
Case 1. p = 2. Using spherical-polar coordinates and the fact that 

/ g{{x, y})doj{y) = <Jd-2 f 9(\xm - t'f''-^'>'Ht, 
Js'i--^ J-i 

where crd_2 is the surface area of S**"^, it follows that 

/•I rd-l(l_^2)M-l/2 rl (1 _ s2)(d-3)/2^^ 

J2 = c / , — / , dr. 

Jo (n-i + Vl - r'^fi^ J-i (1 -|-narccos(rs|a:| + y^l - |a:| Vl - ))2<^ 



Write briefly F{r,t) := 1/[1 + narccos(i|a;| + ^1 - |a;| Vl - r^)]^'^. Then 

(4.20) J2=c ; Il_ / F(r,r.)(l-.^)('^-3)/2^sdr. 

Next, wc apply the substitution u = rs, then switch the order of integration, and 
finally substitute t = Vl — r'^. This gives 

/ J _ / u)(r2_^2-)(d-3)/2^^^^ 

J-iJ\u\ (n-1 + VI - H)2/' 

.1 /.Vl-«^ ^ ^2Mn_^2_ 2N(d-3)/2 

= c/ / f(\/r^,M) ' dtdu. 

J -Jo ' (n-i+i)"^ 

Using the trivial inequality t/{t + n~^) < 1 we conclude that 

1 /.\/T^^ 



J2 < c / / F{\/l-t^, u){l -f- u^)^'^~^^/^dudt. 



1 Jo 



Since ^ sin^/2 ~ ^/l — cos6 for < ^ < tt, we have 

F(\/l-^^^^) ~ {l + n\J I - u\x\ - t^/l - \x\'^^ , < t < Vl - m^- 

But 1 — u\x\ — t\/l — |.'r[2 > if — Vl — <t<Q. Therefore, we can enlarge the 
domain of integration to obtain 



J2 < c 



(1 -e- u^)<^''-^y^dudt 

B2 



(1 + nyT^^u|xP^77r^l^) " 

Here B'^ is the unit disk in K^. We now change the variables {u,t) (a, b), where 



= ^^/T^^Jx^ t + \x\u, b = —\x\t + — |a;|2 u. 
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It is easy to see that this is an orthogonal transformation so that dadb = dudt. 
Hence 

f (1 - a2 _ 52)(rf-3)/2 

J2 < c / , ' — dadb 

- Jb^ (l + nVT^)2'^ 

/•I (l-a2)(d-2)/2 

~ 7-1 (i + 7iyr^)2^ 

< ^ + c / { . da 



„2(T 



/o (1 + nVl - ay 



~ 7l2'^ ^ n<i (1 + s)2<T - „d ' 

since 2(7 > d. Thus (|4.19p is estabUshed when p — 2. 

To prove (|4.19p when p ^ 2 we wiU need the inequaUties 



V2(l + nd(x,2/)) - ^ 

(4.21) < V2{^l^\x\^ + n-')il + nd{x,y)), x,yeB'', 

which foUow readily from (|4.8p . From this and the definition of yV^(a;;n) in (|4.9p 
we get 

(4.22) cW^(n;a;)(l+nd(a;,y))-2A' < W^(n;y) < c>V^(n; a;)(l + n d(a;, y))^^ 
Case 2. < p < 2. Using ((i?^ we obtain 

and hence 



W^(n;y)P/2(l + nd(x,j/))-P - ^^'^ W^(n; y)(l + nd(x, y))- ' 

where r := (ct — 2/i(l/p— l/2))p. By the hypothesis of the lemma t > d. Then the 
above inequality and (|4.19p with p — 2 imply (|4.19p in this case. 

Case 3. 2 < p < oo. Similarly as above by (|4.22p 

VV^,{n,y) - yVf,[n,y)yv^[n,y) ^ (i + y))2^(p/2-i) • 



Consequently, 



W,{y)dy <eW,(n;x)i-^/^ ^ ^'^^^^^^ 



Js^W^{n;y)P/^l + nd{x,y)rP - j W^{n-y){l + nd{x,y)Y ' 

where this time t [a — 2fi{l/2 — \/p))p. Since t > d, the above inequality and 
(|i?Tg|) with p = 2 imply (jiTTg)) in the case 2 < p < oo. □ 

It will be vital for our further development that L!^{x,y) is a Lipl function in 
X (or y) with respect to the distance d{-, •). Throughout the rest of the paper, we 
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denote by B^{r) the closed ball centered at ^ of radius r > with respect to the 
distance c?(-, •) on B'^, i.e. 

B^{r) := {xeB'^ : d{x, ^ <r}, B'^ , r > 0. 

Proposition 4.7. Let £,,y e B'' . Then for all x,z & B^{c*n''^) (c* > 0,n > 1) 
and an arbitrary k, we have 



(4.23) \L^^{x,y)-L^^[U)\<Ck- 



n^+M(x,C) 



VW^(n;y)VW^(n;z)(l + 2))^= ' 

where Ck depends only on k, fi, d, a, and c* . 



Proof. Let ^ > 0. We will use the notation t{x, y; u) := {x, y)+u^l — \x\'^y/l — 
introduced in (HTU). By it follows that 



mx,y)-m,y)\ 



(4.24) 



< c 



< c 



y; u)) - L^{t{C, y; u)) (1 - u^-^du 



dL^J-) \t{x,y;u)-t{^,y;u)\{l-u^r-'du, 



where df — f and /,( is the interval with end points t(x, y\ u) and t(^, y; w). 
As in the proof of Theorem 14.21 by estimate (|4.4p with r = 1 it follows that 



(4.25) 



< Cfcn^^+^ max (1 + n^/l - t] 



using the fact that (1 + n\/\ — t) is an increasing function of r. 



By the definition of t{x,y;u) it follows that (recall x' :— (x, \/l — |a;p)), 
\t{x,y;u) - t(C,y;w)| 



<\W,y')-{^',y')\ + \i-u\Vi^ Vi^-Vi^' 

< I cos d(x, y) ~ cos d(e, y)| + \/2 1 1 - u| ^l-\y\^d(x, 0, 

where we used inequality (|4.8p from Lemma l4. II Denote briefly a := d{x,y) and 
/3:=d(e,y). Then 

I cosd(x,?/) - cosd(^,?/)| = 2 sin ^" ^ ^U in " ^ ^ < i|a - /3|(a + /3) 

< ^\dix,y) - di^,y)\id{x,y) + d{ty)) 

< ^d{x,0{d{x,y) + d{ty)) 



<d{x,^)id{y,z)+c*n-^) 



for z E B^{c*n ). Hence 



|t(x, y; u) - t{^, y;u)\< d{x, 0{d{y, z) + c*n-^) + ^2 |1 - u\^l ~ \y\^d{x, 0- 
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We use this and (I4.25P in ()4.24|) to obtain 



\k 



-1 (1 + ny/1 - t{x,y\uY J -I (1 + n^J\ - t{^,y;u) 

/ {l-u){l-u^y-^du {l-u){l-u^Y-^du 
W-i (l + n^l -t{x,y;u)'' J^i (1 + n^/l -t{^,y;u)'' 



By Lemma 14.41 with 7 = — 1, we have 

Ai<cn^^+^d{x,^){d{y,z)+c*n-~^) , , 

' ^W^{n-x)^W^{n-y){l + nd{x,y)Y 

with a :— fc — 3/i — 1. Note that for y E B'^ and aU z E B^{c*n^^), we have 
l + nd{z,y) ~ l + nd{^,y) and i/l - [zp + c*n"i ^1 - l^p + using g 
Consequently, 

(4.26) ^^'^"^^(^'^^ 



(7-1 



+ nd{y, z)) 

We similarly obtain the same bound for ^2- 

To estimate At, we employ Lemma 14.41 with "f — IJ, and obtain 



A,<cn^^+^^l^d{x,0 

(4.27) 



with (T := fc — 3/i — 4. By cancelling appropriate terms we conclude that ()4.26p 
holds for A3 as well. Exactly in the same way one can see that A4 also satisfies 
(|4.27p and hence (|4.26p . The proof of the proposition is complete. □ 

Operators. We next use the localized polynomials Lt^^ as kernels of linear operators 
defined by 

(4.28) {C^f){x):^b^J f{ymx,y)W^{y)dy, ^ > 0. 

J B-i 

Let E^{f)f^^p denote the best approximation to / £ L^, where := LP{B'^,Wfi), 
from the space HJ^ of all polynomials of degree at most n, that is, 

i?,f(/U :=^^i;nfJl/-g||L.. 

Proposition 4.8. Let a be admissible of type (a). Then the operator satisfies 
the following properties: 

(i) is a polynomial of degree at most 2n; 

(ii) C!^p = p for any polynomial p of degree at most n; 

(iii) for / G LP, 1 < p < 00, 

(4.29) l|/:rjlLS^LS <c and \\C>;j ^ f^.K cE^ if)^,p. 
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Proof. The first two properties are obvious from tlie definition of Since £^ is 
an integral operator, the operator norms and are both 

bounded by 

max / \L'^ix,y)\W^{y)dy. 

Estimate (j4.18p from ProDOsition l4.5l with p = 1 shows that this quantity is bounded 
by a constant independent of n. Then it follows by interpolation that is a 
bounded operator from into for I < p < oo, which yields (|4.29p . □ 

A result of the same nature holds true for more general weight functions of the 
form hf.{x){l — where hK,{x) is some function invariant under a finite 

reflection group, see pp] , 

5. CUBATURE FORMULA ON B'^ 

Cubature formulae on B''' with weights W^{x), /i > 0, which are exact for all 
polynomials of degree n are valuable from many prospectives. Those with positive 
coefficients are preferred for numerical computation and are called positive cubature 
formulae. In the literature, only a handful of positive cubature formulae are known. 
For our purpose of constructing polynomial frames on B'' (see we will need 
positive cubature whose knots are almost equally distributed with respect to the 
distance d[-, ■) introduced in (14. 6p . To the best of our knowledge there are no such 
cubature formulae available up to now. There is a close relation between cubature 
formulae on the unit ball and those on the unit sphere S'^ [IB] . In the following we 
will follow the method used in [10] for constructing cubature formulae on the unit 
sphere. 

One of the difficulties in constructing cubature formulae on S"^' is the lack of 
uniformly distributed points on B"^ . We shall use as a substitute sets of "almost 
equally distributed points" with respect to the distance •) in B^^ which we 
describe in the following. 

Lemma 5.1. For any < e < tt there exists a partition TZ^ of B'^ consisting of 
projections of spherical simplices and a set C B'^ {consisting of their "centers") 
with the properties: 

(i) B"^ = Ui?,e-R, R 0"-^ ihe sets in Tie do not overlap (i?°ni?2 = if Ri ^ R2) . 

(ii) For each R e 7?.£ there is a unique ^ G A'^ such that B(_(c*e) C i? C -B^(e). 
(in) = #7^e < 0**6-". 

Here c* and c** are constants depending only on d. 

Proof. As we already mentioned the distance d{x,y) {x,y € B"^) is the geodesic 
distance between x',y' G 5^. So, we need to subdivide properly 5^. We first 
divide S'^ into 2'^ spherical simplices analogous to the intersections of with the 
octants in R^. Let Oi be the spherical simplex on which all coordinates of ^ € Oi 
are nonnegative and let 

d+l d+1 

Ti:={^t,e, >0, Y.^j^l}, 

j=i i=i 

where {ej} are the standard unit vectors in M'^+^. If w := (1, . . . , 1), then the map 
x(^) := -Tj^r gives an one-to-one correspondence between Oi and 7i. It is readily 
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seen that for ^, C G Ci 

(5.1) ^d{^x)<KO-x{C)\<2Vdd{to. 

Here | • | denotes the Euchdean norm in M''+^ and d{-,-) is the geodesic distance on 

We set M := \2\/de^^~\ and divide the equilateral simplex 7i into M"^ equal 
equilateral subsimplices of side length L = \/2/A'I. We denote by TZl the set of 
all spherical simplices obtained by applying the inverse map to the simplices 
defined above. We similarly define the set of the "centers" of all spherical sim- 
plices in TZl ^PP^ying the inverse map x"^ to the midpoints of the corresponding 
Euclidean simplices. After these preparations, we define TZl ^ of projections 

onto B"^ of all spherical simplices from TZl and we similarly define . 

It is straightforward to show that an equilateral Euclidean simplex with each 
side of length L contains the ball of radius L/ ^2d{d + 1) centered at its midpoint 
and is contained in a ball of radius < L/v^ with the same center. Then (|5.ip yields 
that the corresponding spherical simplex contains the spherical cap centered at its 
center and of radius L / {2d^J 2{d + 1)) and is contained in a spherical cap with the 
same center and radius < V2dL < 2\/d/M < e. This establishes property (ii) of 
Lemma [5. II for the spherical simplices in TZl- Also, we have 4t^X^ = #7?.^ = M'^ < 
{iVde-^y. 

Repeating this procedure with all other initial simplices, we establish the exis- 
tence of the desired partition TZ^ ■ □ 

Definition 5.2. A set X^ C B'^ which, along with an associated partition TZ^ of 
B'^ , has the properties of the sets X^ and TZ^ of Lemma \5.1\ will he called a set of 
almost uniformly e-distributed points on B"^ . 

The following lemma contains additional information about the partition TZ^. 

Lemma 5.3. Let TZ^ be as in Lemma \5.1\ Then for any ^ Xg 

(5.2) \R^\ / Idx^e",/!^ 
and 

(5.3) m,{R^) -.^ w,{x)dx^e\l-\S.\''Y=e''^^^e\^l^ + ef^. 

Here the constants of equivalence depend only on d and fi. 

Proof. To prove ()5.2p we use property (ii) in Lemma 1 5 . II which yields 

(5.4) - |5^(e)| and d{^,dB'^) > c* e. 

We can assume without loss of generality that ^ lies on the positive xi-axis, i.e. 
^ = (Ci,0, ... ,0) and < < 1. The boundary dB^{e) of B^{e) is given by the 
equation a;i^i + -^/l — |a;p-\/l — S^f — cose. A simple manipulation of this identity 
shows that dB^{e) is the ellipsoid 

(xi - ^1 cose)^ 2 2-2 
^ ^_|^|2 +xl + --- + xj = sm^ e. 
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From this it follows that ~ e'^-\/l — (using that sine ~ e) and then 

follows. 

We now turn to the proof of (|5.3p . There are two cases to be considered. 
Case 1. fj,> 1/2. Denote := n {x e B''- : \x\ < |^|}. It is easily seen that 
\R^\^\R^\^e^^T^. Then 

W^ix)dx> [ W,ix)dx>W^{0\R^\--W,{Oe'V^^^^e''{l-\^\^r. 



Since ^ is in the center of TZ(^ by construction, we have -^/l — > ce. Hence, for 
X E R^ C B^{e), inequality ()4.8|) shows that 



which yields 

w,{x)dx < cw.iom - w.iOe'V^^lw = ^'(1 - \er- 



Case 2. < ^ < 1/2. Denote i?+ i?^ n {x e B'^ : |a;| > |^|}. Proceeding as 
above we again get (|5.3p . 

Finally, using (|5.4p we obtain v^l - ICP > sinc*e > ce which implies the last 
equivalence in (|5.3p . The proof of the lemma is complete. □ 



Theorem 5.4. There exists a constant > {depending only on d) such that 
for any n > 1 and a set of almost uniformly e-distrihuted points on B'^ with 
e :— c? jn, there exist positive coefficients \}^{\^ex^^ such that the cubature formula 

[ fix)dx^ ^diO 

is exact for all polynomials of degree < n. In addition, 

with constants of equivalence depending only on /i and d. Here m^{E) := Jg W^{x)dx. 

Note that when — Q the cubature formula of Theorem 15 . 41 can be derived from 
the cubature formula on 5"'+^ from [Tni[IIl[I3] by applying [T51 Theorem 4.2]. 

For the proof Theorem l5.4l we will utilize the idea used in [101 E] (see also [T3]) 
for the construction of a cubature formula on S*"*. 

Assume that (with associated partition TZ^) is a set of almost uniformly e- 
distributed points on B'^ (see Definition 15. 2p . where e = 5/n with n > 1 and & will 
be selected later on. We introduce the following weighted ^i-norm for functions 
defined on B'^: 

(5.5) \\.fhl(x,) E \fiO\mM. 

Also, recall the notation ||/||li = \\f\\L^w^,B'') /gd \fix)\Wf,{x)dx. 
We need a couple of additional results. 

Lemma 5.5. If g E Ilfj, then 



(5.6) ||.g|Ui - ||.9l|£i(^.)| ^ E / \9ix) - gmW^{x)dx < c*5\\g\ 



LI 
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and hence 

(5.7) (1 - c*S)\\9hi^ ^ < (1 + c*S)\\gh^, 

where c* depends only on d and ^. 

Proof. Let £^ be the operator defined in (14.281) . By Proposition 14.51 we have g — 
Ct^g. Using this and the fact that TZ^ is a partition of B'^ (see Lemma [5.ip . we 
obtain 



\\9\\Ll-\\9\\ilix,)\<J2 I \9{x)-g{0\W^,{x)dx 



J2 f f \mx,y)-L'^i^,y)\\giy)\W^iy)dyW^{x)dx 

< 11.911 Li sup ^ / \L>;,{x,y)-L>;,{^,y)\W^{x)dx. 

By Proposition 14. 71 with z = x, it foUows that 
\L>;^{x,y)-L>;^{^,y)\W^{x)dx 

Ckn''+^dix,OW^{x)dx 



< 



Ir^ y/W^{n; x)^W^(n; y)(l + nd(y, x))^ 

Choosing k sufficiently large (fc > d + fi will do) we apply Lemma 14.61 with p = 1 
and use that d{x,^) < Sn~^ for x E to obtain 



sup I \m^,y)-K{U)\w^,{x)d2 



Jb-^ x)^W^{n] y){l + nd{y, x)Y 

The lemma follows. □ 

The Farkas Lemma. A variant of the well known in Optimization Farkas lemma 
will play an important role in the proof of Theorem 15.41 



Proposition 5.6. Let V he a finite dimensional real vector space and denote by 
V* its dual. Let ui,U2,...,Un G V* and suppose u £ V* has the property that 
u(x) > for all x (z V such that Uj{x) > for j = 1, 2, . . . , n. Then there exist 
o-j ^ 0, j = 1, 2, . . . , such that 

n 

(5.8) u = ^ajUj. 

For the proof of this proposition, see e.g. [T]. 

Proof of Theorem \5.4\ First, we choose S := where c* is the constant from 
Lemma [5.51 In applying Proposition 15. 6i we take V := and {%} to be the set 
of all point evaluation functionals {'^^j^eA'^- 

Let the linear functionals u and be defined by 

uig) ■= / g{x)W^{x)dx and u^{g) := u{g) - 7 5(C)"^m(-R?)- 
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Since c*6 = 1/3, the left-hand-side estimate in (|5.7|) yields 

(5.9) llglUi < (3/2)||g||,i(;,,), gent 

Suppose g eU^ and 5(C) > for all ^ & X^. Then using ((5?6| with c*(5 = 1/3 and 
(15.91). we obtain 



= E / lff(^)-5(0|W^M(^)'^a;<c'^^||g|Ui <(l/2)||g||,i(;,^) 

and hence u{g) > (l/2)||g||£i f;^^). Choose 7 := 1/3. Then since g{(,) > 0, ^ G 
we obtain 

u^{g) - - (l/3)||g||,i(;,^) > (l/6)||5||,i > 0. 
Applying Proposition 15.61 to u^, there exist numbers > 0, ^ e Xg, such that 

and hence 

= E («? + (l/3)™^(fi«)).9(0 =: E 3 e 

Therefore, the linear functional X^jgat provides a cubature formula exact 
for all polynomials of degree n. 

Clearly, > m^(_R^)/3 and the estimate < cm^{R^) follows from Lemma [5731 

and Proposition 15 . 71 below. □ 

The last ingredient in bounding A^ from above is the following general result 
that is of independent interest. 

Proposition 5.7. // a positive cubature formula 

(5.10) / .f{x)W^{x)dx ^ E ^«/(^)' ^« > 0' 1^1 < 1' 

is exact for all polynomials of degree < n, then 



(5.11) Xi<cn-''W^{n;0^cn-^{VT^W + n-')^'', U X„ 

where c > depends only on fj, and d. 

Proof RecaU the kernel Kj^iW,,; x, y) defined in (fO)l . Evidently iC„(W^; ^, C) > 
and 

[KUW,,-x, v)f W^(v)dy - K^{W^-x, x). 

B'' 



Let m = [»T./2J. Then it follows that 



KmiW^;tx) 



1 2 



Wfj,{x)dx = 



Hence, the stated result is a consequence of an upper bound for [KmiW^; x, x)] ^ , 
to be established in Proposition 15.91 below. □ 



In order to establish the needed upper bound for [Kn{W^] x, x)] ^ we now con- 
struct a family of well localized polynomials. 
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Lemma 5.8. For any k,in > 1 and ^ E B'^ there exist a polynomial G ^2km 
and a constant c* > depending only on k and d such that P^{^) = 1 and for 
0<-f <k, X e B"^, 



(5.12) 0<P^{x)< — 



< 



(l + mrf(e,x))2'= - (^l-|xP+m-i)7(l+md(5,x))* 



Proof. Let q(0) :~ ^ m"sin(e/2) ) ■ Evidently, q is a trigonometric polynomial of 
degree less than km, g(0) = 1, and 

(5.13) o^^(^)^(rT^' 

For < a < TT, we define the algebraic polynomial Qa{t) by 
n I a\ q{0 ~ a) + q{0 + a) 

TT^) ■ 

It is readily seen that degQa < km, Qa{cosa) — 1, and 

(5.14) 0<Q„(cosg)< I O<0<7r. 

(1 + m\0 — a\Y^ 

Also, Q7r/2 is even and 

(5.15) 0<QV2(i)< (,^^|^^J^,_^/2|)2^ - {l + lw^^ 

Without loss of generality we may assume that ^ = (^i, 0, . . . , 0) with < < 1. 
We choose a € (0,7r/2) so that ~ cos a. Then (|5.14|) gives 

(5.16) »^Q"W^ (l + ^4i,t))2^ ' 

where di(^i,t) := arccos(^it + -^/l — ^^Vl — i^) is the univariate version of the 
distance d{-,-) (see (|4.6p ). We define 

P^{x) Qa{xi)QT,/2(^\J X2 H ^a;^^ 

Clearly, e n^,,„, P^(0 = 1, and by (I5l5l)-(I5IS1) 

(5.17) 0<P^(2;)<-- , — ---rr, xeB'^, 

^ [(l + m|a;,|)(l + mrfi(a,a;i))]2'=' 

where := (a;2, ■ ■ ■ ,Xd) and :== (2:2 + h -t^)^/^. 

It remains to show that P^ obeys (I5.12p . To this end we fist show that 

(5.18) di^,x)<2{\x,\+d,ii,,x,)). 
Denote briefly Xo '■= {xi, 0, . . . , 0). We have 

d{^,x) < d{£,,Xo) + d{xo,x) = di{£_i,xi) + d{xo,x). 
Our next step is to prove the inequality 

(5.19) d{xo,x) <2\x^\. 
Evidently, 

d{x^,x) = arccos((x^, a:')) — arccos (^x^ + \J 1 — x\^ 1 — x\ — x\ — ■ ■ ■ ~ . 
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One easily verifies the inequality arccosi < 2V1 — t, < t < 1, and hence (|5.19p 
will be established if we show that 

1/2 



Denote briefly a := \J\ — x\ and h :— \x.j,\. Then the above inequality is equivalent 
to — a\J 0? — 6^ < }?■ or aV a? — IP- — [c? — b^) > 0. But the latter inequality is 
apparently valid since 



aVa ~b^ - (a - 6 ) = . > 0. 

Thus ([Ogi) is established and hence ([5T^ holds. Combining ([gT7|) with ([OB)) 
gives 

which is the first estimate of P^{x) in (|5.12[) . 

To prove the second estimate in ()5.12p we need the estimate 



5.21 — — < — ^ , xeS'', 

which apparently follows by inequality ()4.8|) in Lemma [4Tl (see also ()4.2ip ). 

Finally, applying (|5.21|) in ()5.20|) we get the second estimate in ()5.12p . which 
completes the proof. □ 

The function A„(a;) [Kn{Wf^;x,x)]^^ is the so called Christoffel function, 
which has the following characteristic property [4, p. 109]: 

(5.22) A„(a;)= min / [Piy)]^Wf,iy)dy, x^B^. 



The localized polynomials in Lemma 15.81 give an upper bound for the Christoffel 
function, used in the proof of Proposition 15.71 

Proposition 5.9. For any /i > and d > 1 there exists a constant c > such that 

(5.23) An{x) <cn-'^W^,{n;x), xeB'^,n>l. 

Proof. Write k := [ma.x{d/2, fi}] + 1 and let n> Ak (the case 1 < n < 4fc is trivial). 
Set m := [n/2k]. By Lemma [5.81 there exists a polynomial Px{y) S HJ^ such that 
Px{x) — 1 and (|5.12p holds with 7 = /i and £,,x replaced by x,y. Then by (|5.22p . 
(|5.12p . and Lemma l4!6l with p = 2, we infer 



An(a;)< / [PAy)rw^{y)dy<c 



Wf,{m;x)Wf,{y)dy 



d W^{m;y){l + md{x,y))'^'' 
< cm^"^Wf^{m]x) < C7i"''yv^(n;a;). 

□ 



For the construction of our frames, we will need the following result which is an 
immediate consequence of Lemma 15.11 and Theorem 15.41 
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Corollary 5.10. There exists a sequence {Xj}fLo of sets of almost uniformly ej- 
distributed points on B'^ [Xj :— X^.) with Ej :— c^2^-'^^ and there exist positive 
coefficients {A^}^g;tj such that the cubature 

(5.24) / fix)W,ix)dx^Y.^if(^'^ 

is exact for all polynomials of degree < 2-'+^. Moreover, ~ m^(B^{2^^)) and 
ifXj ~ 2^'^ with constants of equivalence depending only on d and /i. 

6. Tight polynomial frames (needlets) in L^(_B'',Wp) 

We will apply the general scheme, described in ^ for construction of polynomial 
frames in -L^ := L'^{B'^, W^/j)- To this end we will utilize the localized polynomials 
from Theorem 221 and the cubature formula from Theorem l5.4l (see Corollarv l5.10|) . 

Let a satisfy the conditions: 

(6.1) ae C°°[0,cx)), a>0, suppa C [1/2, 2], 

(6.2) a{t)>c>0, if t e [3/5,5/3], 

(6.3) a'^{t)+a'^{2t) = 1, if [1/2,1]. 

For the construction of such functions, see e.g. [13j. 

We introduce a sequence of polynomial "kernels": Lo{x,y) :— 1 and (see ^ 

CO 

L,{x,y) ■.^J2a(^^y^{W,,;x,y), j = l,2,.... 

i/=0 

We now define the needlets (frame elements) by 

ip^{x) -.^^ ■ Lj{x,^) for ^eXj, j = 0,1,..., 

where Xj is the set of the knots and the A^'s are the coefficients of the cubature 
formula (|5.24p from Corollarv l5.10l We write X := UJ^qXj (see and define the 
needlet system by 

* := {^i}iex- 

Denoting 

(L, * f){x) := / L,{x,y)f{y)W^{y)dy, 
Jb<' 

we get as in (|2.7p the semi-discrete needlet decomposition of L^: 

oo 

f = Y,L,*L,*f for fell. 

The next theorem shows that the needlet system ^I' is a tight frame in L^. 
Theorem 6.1. /// e Ll, then 

OO 

and 
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This theorem foUows at once from Proposition 12 . 21 and Theorem 15.41 (see Corol- 
lary EIOl). 

We finally show that each needlet ■0^ has faster than any polynomial rate of 
decay away from its center (pole) ^ with respect to the distance d(-, •) on B'^, which 
prompted us to coin their name. This property of the needlets is critical for using 
them for decomposition of spaces other than L^. 

Theorem 6.2. For any fc > there exists a constant Cfc > depending only on k, 
II, d, and a such that for ^ £ Xj, j = 0,1, . . . , 

(6.4) \i;,(x)\<Ck , , xeB'^. 

^ ' "^^^ ^' - '^W,i2i;x)il + 2idix,0)'' 

Proof. Estimates (|6.4p follows readily from (|4.10p (see Theorem 14. 2p . taking into 
account that < c2-J'*>V^(2-'; ^ for ^ e X^. □ 

Remark 6.3. Estimate ( |6'.^p and Lemma \4-6\ with p = 2 yield HV'cIIl^ < c for all 
^ £ X, which shows that estimate <\6.4^ is sharp (in a sense). 

7. Appendix 

Proof of Theorem \3.2\ {see We will need the following well known estimate for 
Jacobi polynomials [17, (7.32.6), p. 167]: For a,/? > —1/2 and n>l, 



(7.1) |Fi"'«(cos0)| <c(a,/3) 



mm 
min 



in{n°',n-^/^0-°'-^/^} if < 6* < 7r/2, 

in{n^, n~i/2(7r - 9)-^-^/^} if n/2<e< vr. 



We first prove p.6p for r = 0. We may assume that n > 2k and k > 2, since 
p.6p is trivial when n < 2k (c^ may depend on fc). 

Case 1. < 9 < n/n. Since T{j + a)/T(j + 1) ~ j""-^ and a > /3 > -1/2, it 
follows by ((7T|) that 



|L^'^(cos0)| <c^j2"+i< 



cn^'^+\ 



which yields for < 6* < 7r/n. 

Case 2. vr/n < < vr. A key role here will be played by the identity [T71 (4.5.3), 
p.71]: 

^ (2^ + g + fc + /3 + l)r(^ + a + k + p+\) 



u=0 



_ r(n + o + fc + 1 + /? + 1) (^+k+i.0) 

r(n + d + i) " ^'^>- 



\n + /? + 1) 

Applying summation by parts to the sum in ()3.5p (using ()7.2p with fc = 0), we get 



r(j + Q +1 + + 1) p(a + l,/3). 



r(j + /3 + i) 

We now define the sequence of functions {Ak{t))fLf^ by 



Aoit) := {2t + a + f3 + l)a(^) 
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and inductively 
(7.4) Ak+i{t) := 
It is readily seen that 
(7.5) 



Ak{t) 



Akit + 1) 



2t + a + k + i3 + l 2t + a + k + l3 + 3 



fc > 0. 



and hence supp Afc C [n — fc, 2n] C [n/2, 2n], 1 < fc < n/2. 

Applying summation by parts fc times starting from (13. 5p (using every time 
(|7.2p '). we arrive at the identity: 



(7.6) L^^''{x)=cOJ2A,{j) 



T{j + a + k + (3+ 1) p{a+k,f3) j.^^ 



By (|7.5p it readily follows that 

(7.7) II^^^IU f^^-^-'lla^^+^'llcc. ™>0 

and inductively, using (|7.4p . it follows that 



(7.8) 



\A 



(™)|i 

k Woo 



< cn-™-''=+i max Ha^'^Hoo, m > 0, fc > 2, 



0<u<m+k 



where c = c(fc, m, a, (3). 

Now, from ([7?T|) and ((7?5)) with m = 0, we obtain for Ti/n <9< tt/2 



2n 



|L^"3(COS0)| <C ^ ^-2fc+l^-a+fc-l/2^-a-fc-l/2 



j=n-k 



t2q+2 



< cn2"+2(n^^)-"-^-l/2 < c-i^ — , 



and for 7r/2 < 6* < tt 



-fe+2Q+2 



,2a+2 



< C 



(1 +7l6 



|L^''^(cos6')| < c ^ n-2fe+ij-"+fe/ < 

Thus p.6p is established when r ~Q. 

The case when r > 1 is an easy consequence of Markov's inequality: If Q € n„i, 
then liQ'lU-M] < 2m2(6-a)-i||Q|Uoo[,,6]. 

We give the proof for r = 1 only; in general it follows inductively. Clearly, p.6p 
with r = is equivalent to 

(7.9) \L':^P{x)\<c—^^^^-, xehl,!]. 

(1 + rtVl — 

If X e [0, 1], then by (17. 9p and Markov's inequality 



< 



< 



L=[-l,: 
2a+4 



8n2(l + a;)-i||i^^^|Uoe[_i,,] 



< c- 



(1 + n^/T^)'= 

which is (|3.6p with r = 1. For x £ [—1,0) we apply Markov's inequality on [—1, 0] 
which leads readily to the same result. The proof of Theorem 13.21 is complete. □ 
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